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ABSTRACT 

A review is made of recent magnetohydrodynamic (MHD) theory and simulations of origin of jets 
from accretion disks. Many compact astrophysical objects emit powerful, highly-coUimated, oppositely 
directed jets. Included are the extra galactic radio jets of active galaxies and quasars, and old compact 
stars in binaries, and emission line jets in young stellar objects. It is widely thought that these different 
jets arise from rotating, conducting accretion disks threaded by an ordered magnetic field. The twisting 
of the B field by the rotation of the disk drives the jets by magnetically extracting matter, angular 
momentum, and energy from the accretion disk. Two main regimes have been discussed theoretically, 
hydromagnetic winds which have a significant mass flux, and Poynting flux jets where the mass flux is 
negligible. Over the past several years, exciting new developments on models of jets have come from 
progress in MHD simulations which now allow the study of the origin - the acceleration and coUimation 
- of jets from accretion disks. Simulation studies in the hydromagnetic wind regime indicate that the 
outflows are accelerated close to their region of origin whereas the coUimation occurs at much larger 
distances. 

Subject headings: jets, accretion disks — outflows: jets — galaxies: magnetic fields — plasmas — stars 



1. INTRODUCTION 

Powerful, highly-coUimated, oppositely directed jets are 
observed in active galaxies and quasars (see for exam- 
ple Bridle & Eilek 1984), old compact stars in binaries 
(Mirabel & Rodriguez 1994), and emission line jets in 
young stellar objects (Mundt 1985; Biihrke, Mundt, & Ray 
1988). A broad spectrum of ideas and models have been 
put forward to explain astrophysical jets (see reviews by 
Begelman, Blandford, & Rees 1984 and Bisnovatyi-Kogan 
1993). The matter is thought to go to the jet from the 
inner region of an accretion disk surrounding the compact 
object - a star or black hole. The disk matter must then 
be accelerated to a velocity higher that the escape veloc- 
ity from the central object. Further, the jet matter should 
have sufficient momentum to propagate through surround- 
ing inter-stellar or intra-galactic matter out to large dis- 
tances. 

An ordered magnetic field is widely thought to have an 
essential role in jet formation from a rotating accretion 
disk. Two main regimes have been considered in theo- 
retical models, the Poynting flux regime where the energy 
outflow from the disk is carried mainly by the electromag- 
netic field with the energy carried by the matter small, 
and the hydromagnetic regime where the energy is carried 
by both the electromagnetic field and the kinetic flux of 
matter. Poynting flux models for the origin of jets were 
proposed by Lovelace (1976) and Blandford (1976) and 
studied further by Lovelace et al. (1987) and Colgate & 
Li (1998). In these models the rotation of a Keplerian ac- 
cretion disk twists a poloidal field threading the disk, and 
this results in outflows out of the disk which carry angular 



momentum (in the twist of the field) and energy (in the 
Poynting flux) away from the disk, thereby facilitating the 
accretion of matter. 

Important questions to be answered by the jet models 
include the following: 1. What is the main driving force 
pushing matter into the jet? 2. What determines the mass 
flow rate in the jet Mj, and what fraction is this of the ac- 
cretion rate? 3. What physics determines the asymptotic 
speed or Lorentz factor of the bulk flow? 4. What de- 
termines the coUimation of the jet and at what distance 
from the central object does the jet become collimated? 5. 
What is the acceleration mechanism of leptons to Lorentz 
factors 7 ^ 10^ — 10"^ in the radio jets? Observations of 
Blazars indicate that 7 ^ 10"^ — 10^ in some objects. 

The focus of recent work has been the hydromagnetic 
regime of jet formation for the geometry sketched in Figure 
1. A strong case for hydromagnetic jets as an explanation 
of jets in protostellar systems emerges because the temper- 
ature of the inner regions of these systems is insufhcient to 
permit driving by thermal or radiation pressure (Konigl & 
Ruden 1993). Part of the investigations have been analyt- 
ical or semi-analytical and outgrowths of the self-similar 
solution of Blandford & Payne (1982) (Pudritz & Norman 
1986; Konigl 1989; PcUetier & Pudritz 1992; Contopoulos 
& Lovelace 1994). The outflows in this model are often 
referred to as "centrifugally driven" owing to the driving 
force close to the disk: If the poloidal magnetic field lines 
diverge from the disk surface (making an angle with the 
z-axis of more than 30°), then the sum of the gravitational 
and centrifugal forces is in the +z direction for an MHD 
fluid particle which initially tends to maintain a constant 
angular rotation rate. The self-similar models are unsat- 
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Fig. 1. — Sketch of an accretion disk threaded by a magnetic field for conditions which may lead to hydromagnetic jet 
formation. 



isfactory in the respect that they must be cutoff at small 
cylindrical radii, r < r„ii„. This is the most important 
region of the jet flow. Observations of optical stellar jets 
(Mundt 1985) reveal jet velocities ~ 200-400 km/s, which 
are comparable to the Keplerian disk velocity close to the 
star's surface. This suggests that the jets originate from 
the inner region of the disk close to the star (Shu et al. 
1988; Pringle 1989). 

The limitations of analytical models has motivated ef- 
forts to study jet formation using MHD simulations. Sim- 
ulation studies of hydromagnetic jet formation have ad- 
dressed two main regions: The jet formation region where 
the matter enters with sub slow-magnctosonic speed and 
exits with super fast-magnetosonic speed. The second re- 
gion includes the disk and the problem of the Velikhov 
(1959) - Chandrasckhar (1981) - Balbus-Hawley (1998) in- 
stability and the resulting 3D MHD turbulence. A num- 
ber of studies have addressed the coupled problem of the 
disk and near jet regions (Uchida & Shibata 1985; Shibata 
& Uchida 1986; Stone & Norman 1994; BeU and Lucek 
1995;). MHD simulations of the near jet region have been 
carried out by several groups (Bell 1994; Ustyugova et al. 
1995, 1998; Koldoba et al. 1995; Romanova et al. 1997, 
1998; Meier et al. 1997; Ouycd & Pudritz 1997). 

Here, we first review recent results on MHD simula- 
tions of hydromagnetic jet formation and later discuss the 
Foynting flux regime. Section 2 discusses general con- 
siderations of hydromagnetic outflows. Sections 3 dis- 
cusses MHD simulations which give non-stationary and 
stationary hydromagnetic outflows. Sections 4 describes 
the Foynting flux regime which remains to be fully ex- 
plored by simulations. Section 5 gives the conclusions. 



2. BASICS OF MHD OUTFLOWS FROM DISKS 

The main forces which drive a hydromagnetic outflow 
from a disk threaded by a magnetic field are the centrifu- 



gal force and the magnetic pressure gradient force. If disk 
has a hot corona, the pressure gradient may also be impor- 
tant. We neglect the radiative force but in this regard see 
Fhinney (1987). Accreting matter of the disk carries mag- 
netic field inward thus generating a component of the 
magnetic field as sketched in Figure 1. On the other hand, 
rotation of the disk acts to generate a toroidal component 
of the field (< if 5^ > 0). 

For a sufficiently inclined magnetic field {6 in Figure 1 
sufficiently large) , outflows can result from the centrifugal 
force (Blandford & Fayne 1982) and/or the magnetic pres- 
sure gradient force (— B'^/{8n)) (Lovelace et al. 1989, 
1991; Koupelis fc Van Horn 1989). This depends on the 
ratio of energy densities at the base of the outflow at the 
inner radius of the disk denoted r^. Thus, the main pa- 
rameters are eth = {cs/vk)^ and Eb = {vA/vK)f, where 
vk is Keplerian velocity, Cg the sound speed, va the Alfven 
speed, and the i subscript indicates evaluation on the sur- 
face of the disk at its inner radius, r = ri. For £b ~ 1 the 
outflow is magnetically driven, whereas for Sth ^ 1, the 
flow is thermally driven. 

Frocesses in the disk are of course coupled to the out- 
flows (Lovelace et al. 1994, 1997; Falcke, Malkan, k. Bier- 
mann 1995). However, it is difficult to simulate both re- 
gions simultaneously because the time scales of the accre- 
tion and outflow are in general very different. The accre- 
tion is much slower. On the other hand the processes in the 
disk may involve the small scale MHD instability of Chan- 
drasekhar, Velikhov, Balbus, and Hawley, and therefore 
require high spatial resolution. Stone & Norman (1994) 
attempted to simultaneously simulate the internal MHD 
dynamics of a disk and MHD dynamics of outflows. This 
proved impractical because essentially all of the spatial 
resolution was needed for treating the unstable dynamics 
of the disk. Also, there was the problem that the initial 
configuration was far from equilibium. Simulation of the 
internal MHD disk dynamics has led several groups to the 
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Fig. 2. — Simulations results for stationary MHD outflow obtained using a spherical coordinate system (Ustyugova et al. 
1998). The solid lines represent the poloidal magnetic field, and the arrows the velocity vectors. The dashed lines in the 
left-hand plot represent the slow magnetosonic surface (the lowest dashed line), the Alfven surface (the middle line), and 
the fast magnetosonic surface (the top line). In the right-hand panel, the dashed lines are surfaces of constant toroidal 
current density, while the lines on the outer boundary are the projections of the fast magnetosonic Mach cones. 



problem of simulating 3D MHD turbulence in a sheared 
flow of a patch of a disk (for example, Hawley et al. 1995; 
Brandenburg et al. 1995). This is a much larger project 
than that of understanding MHD outflows. At the same 
time it is widely thought, and observations of cataclysmic 
variables support the view, that the disk turbulence - in- 
cluding MHD turbulence - can be modeled approximately 
using the Shakura (1973), Shakura & Sunyaev (1973) "al- 
pha" viscosity model (Eardlcy & Lightman 1975; Coroniti 
1981). In contrast with the internal disk dynamics, there 
is theoretical and simulation evidence that the outflows 
can be treated using axisymmetric (2D) MHD (Bland- 
ford & Payne 1982; Lovelace et al. 1991; Ustyugova et 
al. 1995). Here, we consider outflows from a disk repre- 
sented as a boundary condition. This approach has subse- 
quently been adopted by other groups (Meier et al. 1997; 
Ouyed & Pudritz 1997). This treatment of the disk is jus- 
tified for outflows from a disk where the accretion speed 
is small compared with the Keplerian speed (Ustyugova et 
al. 1995). 



3. NUMERICAL SIMULATIONS OF MHD OUTFLOWS 

In order to test the analytical models of stationary out- 
flows, MHD simulations of flows from a disk treated as a 
boundary condition have been carried out by a number of 
groups. 

3.1. Non-Stationary Outflows 

The initial magnetic field configuration was chosen so 
that the magnetic field was significantly inclined to the 
disk {9 > 30°) over most of the disk surface (Ustyugova 
et al. 1995; Koldoba et al. 1995). The simulations in- 
volve solving the complete system of ideal non-relativistic 



MHD equations using a Godunov-type code assuming ax- 
isymmetry but taking into account all three components 
of velocity and magnetic field. Matter of the corona was 
initially in thermal equilibrium with the gravitating cen- 
ter. At f = 0, the disk is set into rotation with Keple- 
rian velocity and at the same time matter is pushed from 
the disk with a small poloidal velocity equal to a frac- 
tion of the slow magnetosonic velocity ( Vp = aVsm, with 
a = 0.1 — 0.9). A relatively high temperature and small 
magnetic field was considered. We found that at the max- 
imum of the outfiow, matter is accelerated to speeds in 
excess of the escape speed and in excess of the fast magne- 
tosonic speed within the simulation region (~ lOOrj). The 
acceleration is due to both thermal and magnetic pressure 
gradients. The outfiow coUimates within the simulation 
region due to strong amplification, 'wrapping up' of the 
toroidal magnetic field and the associated pinching force. 

However, the outflows arc not stationary. The matter 
flux grows to a peak and then decreases to relatively small 
values. The strong collimation of the outflow reduces the 
divergence of the field away from the z— axis {9 < 30°) 
and this "turns off" the outflow of matter and leads to 
flow velocities less than the escape speed. Thus, this sim- 
ulation is an example of a temporary outburst of matter 
to a jet. Unfortunately, this type of flow has a significant 
dependence on the initial conditions. 

3.2. Stationary Outflows 

More recently, stationary magnetohydrodynamic out- 
flows from a rotating accretion disk have been obtained by 
time-dependent axisymmetric simulations by Romanova et 
al. (1997) and systematically analyzed by Ustyugova et 
al. (1998). The initial magnetic field in the latter work 
was taken to be a split-monopole poloidal field configu- 
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Fig. 3. — Dependences of different velocities on distance s measured in units of from the disk along the second magnetic 
field line away from the axis in Figure 2 (Ustyugova et al. 1998). This field line "starts" from the disk at r « Qn where it 
has an angle 9 « 28° relative to the z— axis. The velocities are measured in units of {GM /riY/"^ . Here, Vp is the poloidal 
velocity along the field line and v±_ is the poloidal velocity perpendicular to the field line. Also, VAp is the poloidal Alfven 
velocity, c/^ is the fast magnetosonic velocity, and Vesc is the local escape velocity. 



ration (Sakurai 1987) frozen into the disk. The disk was 

treated as a perfectly conducting, time-independent den- 
sity boundary [p(r)] in Keplerian rotation which is differ- 
ent from our earlier specification of a small velocity out- 
flow (§3.1, Ustyugova et al. 1995). The outflow veloc- 
ity from the disk is determined self-consistently from the 
MHD equations. The temperature of the matter outflow- 
ing from the disk is small in the region where the magnetic 
field is inclined away from the symmetry axis (c^ -C v|-), 

but relatively high (c^ ^ w|-) at very small radii in the 
disk where the magnetic field is not inclined away from 
the axis. We have found a large class of stationary MHD 
winds. Within the simulation region, the outflow acceler- 
ates from thermal velocity (~ Cg) to a much larger asymp- 
totic poloidal flow velocity of the order of O.h^GM/n, 
where M is the mass of the central object and is the 
inner radius of the disk. This asymptotic velocity is much 
larger than the local escape speed and is larger than fast 
magnetosonic speed by a factor of ~ 1.75. The acceleration 
distance for the outflow, over which the flow accelerates 
from to, say, 90% of the asymptotic speed, occurs at a 
flow distance of about SOrj. The flows are approximately 
spherical outflows, with only small coUimation within the 
simulation region. The coUimation distance over which the 
flow becomes coUimatcd (with divergence less than, say, 
10°) is much larger than the size of our simulation region. 
Close to the disk the outflow is driven by the centrifugal 
force while at all larger distances the flow is driven by the 
magnetic force which is proportional to — V(ri?0)^, where 
is the toroidal field. 



The stationary MHD flow solutions allow us (1) to com- 
pare the results with MHD theory of stationary flows, (2) 
to investigate the influence of different outer boundary 
conditions on the flows, and (3) to investigate the influ- 
ence of the shape of the simulation region on the flows. 
The ideal MHD integrals of motion (constants on flux sur- 
faces discussed by Lovelace et al. 1986) were calculated 
along magnetic field lines and were shown to be constants 
with accuracy 5 — 15%. Other characteristics of the nu- 
merical solutions were compared with the theory, including 
conditions at the Alfven surface. 

Different outer boundary conditions on the toroidal 
component of the magnetic field can significantly infiuence 
the calculated fiows. The commonly used "free" boundary 
condition on the toroidal field leads to artificial magnetic 
forces on the outer boundaries, which can give spurious 
coUimation of the fiows. New outer boundary conditions 
which do not give artificial forces have been proposed and 
investigated by Ustyugova et al. (1998). 

The simulated fiows may also depend on the shape of the 
simulation region. Namely, if the simulation region is elon- 
gated in the z— direction, then Mach cones on the outer 
cylindrical boundaries may be partially directed inside the 
simulation region. Because of this, the boundary can have 
an artificial infiuence on the calculated flow. This effect 
is reduced if the computational region is approximately 
square or if it is sphcaic-al as in Figure 2. Simulations 
of MHD outflows with an elongated computational region 
can lead to artificial coUimation of the flow. 

Recent simulation studies have treated MHD outflows 
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from disks with more general initial B field configura- 
tions, for example, that where the poloidal field has differ- 
ent polarities as a function of radius (Hayashi, Shibata, & 
Matsumoto 1996; Goodson, Winglee, & Bohm 1997; Ro- 
manova et al. 1998). The differential rotation of the foot- 
points of B field loops at different radii on the disk surface 
causes twisting of the coronal magnetic field, an increase in 
the coronal magnetic energy, and an opening of the loops 
in the region where the magnetic pressure is larger than 
the matter pressure (/3~1) (Romanova et al. 1998). In the 

region where f3^1, the loops may be only partially opened. 
Current layers form in the narrow regions which separate 
oppositely directed magnetic field. Reconnection occurs 
in these layers as a result of the small numerical magnetic 



diffusivity. In contrast with the case of the solar corona, 
there can be a steady outflow of energy and matter from 
the disk surface. We find that the power output mainly 
in the form of a Poynting flux. Opening of magnetic field 
loops and subsequent closing can give reconnection events 
which may be responsible for X-ray flares in disks around 
both stellar mass objects and massive black holes (Hayashi 
et al. 1996; Goodson, Winglee, & Bohm 1997; Romanova 
et al. 1998). 

4. POYNTING FLUX JETS 

In a very different regime from the hydromagnetic flows 
discussed in §3, a Poynting flux jet transports energy and 
angular momentum mainly by the electromagnetic fields 




Fig. 4. — Field configuration of a Poynting flux jet from Colgate and Li (1998). 
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with only small contributions of the matter (Lovelace 
ct al. 1987; Colgate and Li 1998). A steady Poynt- 
ing jet - sketched in Figure 4 - can be characterized in 
the lab frame by its asymptotic (large distance) magnetic 
field = —Bo[ro/rj{z)] at the jet's edge, r = rj{z), 
where rg is the jet's radius at z = and Bq is the 
poloidal field strength at this location. The electric field 
in the jet is E = — v x B/c and consequently the en- 
ergy flux or luminosity of the jet is Lj = vB^Tq/S ~ 
4 X lo45crg/s(i;/c)(ro/10i''cm)2(Bo/10^G)2. Propagating 
disturbances in such field dominated jets provide a simple, 
but self-consistent physical model for the gamma ray flares 
observed in Blazars (Romanova & Lovelace 1997; Levin- 
son 1998; Romanova these proceedings). Owing to pair 
production close to the black hole, the main constituent of 
a Foynting flux jet may be electron-positron pairs. 

5. CONCLUSIONS 

MHD simulations carried out by a number of groups 
over the last several years support the idea that an or- 
dered magnetic field of an accretion disk can give powerful 
outflows of matter, energy, and angular momentum. The 
studies so far have been in the hydromagnetic regime and 
find asymptotic flow speeds of the order of the maximum 
Keplerian velocity of the disk. In contrast, observed VLBI 
jets in quasars and active galaxies point to bulk Lorentz 
factors of order 10 - much larger than the disk Lorentz 
factor. This may be a result of the relativistic dynamics 
not included here, but more likely it reflects the fact that 
these jets are in the Poynting flux regime. Also, these jets 
may involve energy extraction from a rotating black hole 
(Blandford and Znajek 1977; Livio et al. 1998). 
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grant AST-9320068 and NASA grant NAG5 631 1'. The 
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96-02-17113. Also, the research described here was made 
possible in part by Grant No. RPl-173 of the U.S. Civilian 
R&D Foundation for the Independent States of the Former 
Soviet Union. 

Discussion 

Duccio Macchetto: What is the coUimation distance in the 
case of relativistic jets? Does it also scale and collimate at 
distances ^ lOOri? If so it would still be too close to be 

resolved spatially. 

Richard Lovelace: Simulations of the origin of relativistic 
jets remains to be done systematically. In the hydromag- 
netic regime, the flows are likely to be similar to those 
in the non-relativistic limit with little collimation at dis- 
tances < lOOrj. In contrast I believe that the Poynting 
flux jets collimate rapidly over distances of order several 
fi- 

Julian Krolik: What parameters determine the mass out- 
flow rate? 

Richard Lovelace: For the hydromagnetic outflows we flnd 
in the non-relativistic limit that roughly Mj oc rfBf, 
where is the inner radius of the disk [(2 — 6)GM/c'^] 
and Bi is the poloidal field at r,. However, Mj will be 
less than some fraction of the disk accretion rate. For the 
Poynting flux jets, Mj is negligible compared with the disk 
accretion rate. 

Michael Dopita: In the wind dominated case, can the wind 
itself become of sufficiently high column density to become 
optically thick to Thompson scattering and thereby form 
a photosphere in the wind? 

Richard Lovelace: A hydromagnetic outflow with Mj say 
one-tenth of the disk accretion rate would have a photo- 
sphere well above the disk. 
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